Using Protein Homology Models for Structure-Based Studies: Approaches to Model Refinement

Homology modeling is a computational methodology to assign a 3-D structure to a target protein when experimental data are not available. The methodology uses another protein with a known structure that shares some sequence identity with the target as a template. The crudest approach is to thread the target protein backbone atoms over the backbone atoms of the template protein, but necessary refinement methods are needed to produce realistic models. In this mini-review anchored within the scope of drug design, we show the validity of using homology models of proteins in the discovery of binders for potential therapeutic targets. We also report several different approaches to homology model refinement, going from very simple to the most elaborate. Results show that refinement approaches are system dependent and that more elaborate methodologies do not always correlate with better performances from built homology models

The statistical-thermodynamic basis for computation of binding affinities: a critical review

Although the statistical thermodynamics of noncovalent binding has been considered in a number of theoretical papers, few methods of computing binding affinities are derived explicitly from this underlying theory. This has contributed to uncertainty and controversy in certain areas. This article therefore reviews and extends the connections of some important computational methods with the underlying statistical thermodynamics. A derivation of the standard free energy of binding forms the basis of this review. This derivation should be useful in formulating novel computational methods for predicting binding affinities. It also permits several important points to be established. For example, it is found that the double-annihilation method of computing binding energy does not yield the standard free energy of binding, but can be modified to yield this quantity. The derivation also makes it possible to define clearly the changes in translational, rotational, configurational, and solvent entropy upon binding. It is argued that molecular mass has a negligible effect upon the standard free energy of binding for biomolecular systems, and that the cratic entropy defined by Gurney is not a useful concept. In addition, the use of continuum models of the solvent in binding calculations is reviewed, and a formalism is presented for incorporating a limited number of solvent molecules explicitly

Effect of HIV infection and antiretroviral therapy on hepatitis B virus (HBV)-specific T cell responses in patients who have resolved HBV infection

Coinfection with hepatitis B virus (HBV) is a common occurrence in human immunodeficiency virus (HIV)–positive patients and an increasing cause of morbidity and mortality. The CD8+ T cell response is critical for long-term control of HBV in patients resolving acute infection. Here, we examine the effect of HIV on HBV-specific CD8+ T cell responses in patients who have resolved HBV infection. A cross-sectional study showed a reduction in HBV-specific CD8+ T cell responses in HIV-positive, HBV-immune patients, compared with those in HIV-negative, HBV-immune patients. A longitudinal study of a subgroup of patients examined whether this attrition could be reversed by effective antiretroviral therapy. The introduction of highly active antiretroviral therapy (HAART) resulted in reconstitution of some HBV-specific CD4+ and CD8+ T cell responses, in association with restoration of CD4+ T cell counts. These data provide a mechanism to account for the observed impairment of control of HBV infection in the setting of HIV infection and support the ability of HAART to reconstitute functionally active T cell responses

Block bond-order potential as a convergent moments-based method

The theory of a novel bond-order potential, which is based on the block Lanczos algorithm, is presented within an orthogonal tight-binding representation. The block scheme handles automatically the very different character of sigma and pi bonds by introducing block elements, which produces rapid convergence of the energies and forces within insulators, semiconductors, metals, and molecules. The method gives the first convergent results for vacancies in semiconductors using a moments-based method with a low number of moments. Our use of the Lanczos basis simplifies the calculations of the band energy and forces, which allows the application of the method to the molecular dynamics simulations of large systems. As an illustration of this convergent O(N) method we apply the block bond-order potential to the large scale simulation of the deformation of a carbon nanotube.Comment: revtex, 43 pages, 11 figures, submitted to Phys. Rev.

Effective interaction between helical bio-molecules

The effective interaction between two parallel strands of helical bio-molecules, such as deoxyribose nucleic acids (DNA), is calculated using computer simulations of the "primitive" model of electrolytes. In particular we study a simple model for B-DNA incorporating explicitly its charge pattern as a double-helix structure. The effective force and the effective torque exerted onto the molecules depend on the central distance and on the relative orientation. The contributions of nonlinear screening by monovalent counterions to these forces and torques are analyzed and calculated for different salt concentrations. As a result, we find that the sign of the force depends sensitively on the relative orientation. For intermolecular distances smaller than $6\AA$ it can be both attractive and repulsive. Furthermore we report a nonmonotonic behaviour of the effective force for increasing salt concentration. Both features cannot be described within linear screening theories. For large distances, on the other hand, the results agree with linear screening theories provided the charge of the bio-molecules is suitably renormalized.Comment: 18 pages, 18 figures included in text, 100 bibliog

Age-related craniofacial differences based on spatio-temporal face image atlases

A number of studies have been developed recently in order to explore associations between craniofacial differences and genetics. Most of these works have been based on spatial face image models, adjusted for the counter effects of age. This approach provides a limited understanding of normal and abnormal craniofacial development owing to the lack of age progression information. Here, the authors propose and implement an imaging framework that combines facial landmark positioning, non-rigid registration, novel age-dependent face modelling and common distance metrics to disclose the most facial differences that vary across the time due to the subjects' age. All the experiments carried out and corresponding results presented here are based on a database comprising ordinary two-dimensional (2D) frontal face images of Down Syndrome (DS) and control sample groups. A number of craniofacial metrics have been successfully identified that highlight statistically significant and clinically relevant differences between the controls and the faces associated with DS within the age range from 1 to 18 years old, producing realistic unbiased face models with similar level of detail at all age-intervals, despite the small sample size available

Adsorption of mono- and multivalent cat- and anions on DNA molecules

Adsorption of monovalent and multivalent cat- and anions on a deoxyribose nucleic acid (DNA) molecule from a salt solution is investigated by computer simulation. The ions are modelled as charged hard spheres, the DNA molecule as a point charge pattern following the double-helical phosphate strands. The geometrical shape of the DNA molecules is modelled on different levels ranging from a simple cylindrical shape to structured models which include the major and minor grooves between the phosphate strands. The densities of the ions adsorbed on the phosphate strands, in the major and in the minor grooves are calculated. First, we find that the adsorption pattern on the DNA surface depends strongly on its geometrical shape: counterions adsorb preferentially along the phosphate strands for a cylindrical model shape, but in the minor groove for a geometrically structured model. Second, we find that an addition of monovalent salt ions results in an increase of the charge density in the minor groove while the total charge density of ions adsorbed in the major groove stays unchanged. The adsorbed ion densities are highly structured along the minor groove while they are almost smeared along the major groove. Furthermore, for a fixed amount of added salt, the major groove cationic charge is independent on the counterion valency. For increasing salt concentration the major groove is neutralized while the total charge adsorbed in the minor groove is constant. DNA overcharging is detected for multivalent salt. Simulations for a larger ion radii, which mimic the effect of the ion hydration, indicate an increased adsorbtion of cations in the major groove.Comment: 34 pages with 14 figure

Novel Druggable Hot Spots in Avian Influenza Neuraminidase H5N1 Revealed by Computational Solvent Mapping of a Reduced and Representative Receptor Ensemble

The influenza virus subtype H5N1 has raised concerns of a possible human pandemic threat because of its high virulence and mutation rate. Although several approved anti-influenza drugs effectively target the neuraminidase, some strains have already acquired resistance to the currently available anti-influenza drugs. In this study, we present the synergistic application of extended explicit solvent molecular dynamics (MD) and computational solvent mapping (CS-Map) to identify putative ‘hot spots’ within flexible binding regions of N1 neuraminidase. Using representative conformations of the N1 binding region extracted from a clustering analysis of four concatenated 40-ns MD simulations, CS-Map was utilized to assess the ability of small, solvent-sized molecules to bind within close proximity to the sialic acid binding region. Mapping analyses of the dominant MD conformations reveal the presence of additional hot spot regions in the 150- and 430-loop regions. Our hot spot analysis provides further support for the feasibility of developing high-affinity inhibitors capable of binding these regions, which appear to be unique to the N1 strain